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Abstract. We prove that under n 1 ' 3 scaling, the limiting distribution as n — > oo of the free energy 
of Seppalainen's log-Gamma discrete directed polymer is GUE Tracy- Widom. The main technical 
innovation we provide is a general identity between a class of n-fold contour integrals and a class of 
Fredholm determinants. Applying this identity to the integral formula proved in [11] for the Laplace 
transform of the log-Gamma polymer partition function, we arrive at a Fredholm determinant which 
lends itself to asymptotic analysis (and thus yields the free energy limit theorem). The Fredholm 
determinant was anticipated in [7] via the formalism of Macdonald processes yet its rigorous proof 
was so far lacking because of the nontriviality of certain decay estimates required by that approach. 



1. Introduction and main results 

The log-Gamma polymer was introduced and studied by Seppalainen [15]. 

Definition 1.1. Let 9 be a positive real. A random variable X has inverse-Gamma distribution 
with parameter 9 > if it is supported on the positive reals where it has distribution 

F(X G dx) = — i-roT 9 " 1 exp j — 1 dx. 

T(9) { x) 

We abbreviate this X ~ Y~ 1 {9). 

Definition 1.2. The log-Gamma polymer partition function with parameter 7 > is given by 

Z(n,N)= Ii di d 

7r:(l,l)->-(n,jV) (i,i)G7r 

where ir is an up/right directed lattice path from the Euclidean point (1,1) to (n,N) and where 

In [15] it was proved that 

logZ(n,n) - varlogZ(n,n) 
lim = j 7 , Inn sup — rz < C 

n.-»oo n rwoo n 2 ' 6 

where / 7 = —2^(7/2) and C is a large constant. Here = [logr]'(x) is the digamma function. 
The scale of the variance upper-bound is believed to be tight, since directed polymers at positive 
temperature should have KPZ universality class scalings (see e.g. the review [10]). Moreover, it is 
believed that, when centered by n/ 7 and scaled by n 1 / 3 , the distribution of the free energy log Z(n, n) 
should limit to the GUE Tracy- Widom distribution [16]. 

We presently prove this form of KPZ universality for the log-Gamma polymer for 7 < 7* for some 
7* > 0. This assumption is purely technical and comes from the asymptotic analysis. It is likely 
that this assumption can be removed following the approach of [8] , where a similar assumption was 
removed in the case of the semi-discrete polymer. For this model 7 plays a role akin to temperature. 



1 
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Theorem 1. There exists 7* > such that the log- Gamma polymer free energy with parameter 
7 £ (0,7*) has limiting fluctuation distribution given by 

m ( log Z(n, n) — n f~ 
lim P 6 v ^ < r 

where / 7 = — 2^(7/2), g 7 = — 2 x &"(7/ / 2) and -Fgue Gf/E 1 Tracy-Widom distribution function. 

We give the proof of this theorem in Section 2. There are two ingredients in the proof, and 
then some asymptotic analysis. The first ingredient is the ra-fold integral formula given in [11] for 
the Laplace transform of the polymer partition function. This is given below as Proposition 1.4. 
The second ingredient in the proof is a general identity between a class of n-fold contour integrals 
and a class of Fredholm determinants. This is given below as Theorem 2. Applying this identity to 
Proposition 1.4 yields Corollary 1.8 which is a new Fredholm determinant expression for the Laplace 
transform of the log-Gamma polymer partition function. This formula lends itself to straightforward 
asymptotic analysis, as is done in Section 2. 

The log-Gamma polymer may be generalized, as done in [11], so that the distributions of the 7$ 
depend on two collections of parameters. 

Discrete directed polymer partition functions, under intermediate disorder scaling [2, 13], converge 
to the solution of the multiplicative stochastic heat equation (whose logarithm is the KPZ equation). 
If the two collections of parameters determining the distributions of the 7jj are tuned correctly, 
then the initial data for the limiting stochastic heat equation is determined by two collections 
of parameters as well. The Fredholm determinant formula of Corollary 1.8 should limit to an 
analogous formula for the Laplace transform of the stochastic heat equation with this general class 
of initial data which would be a finite temperature analog of the results of [9]. When only one 
of the collections of parameters is tuned, this formula was computed in [8] via a similar limit of 
the Fredholm determinant formula for the Laplace transform of the semi-discrete polymer partition 
function (see also [12] in the case where only a single parameter is tuned), and this is a finite 
temperature analog of the results of [6]. 

Definition 1.3. The Sklyanin measure sn on is given by 

1*1* 

(i- 1 ) s N{dw U . . . , dWN ) = n r( } n 

v ; i,j=l v * ■>' i=l 

The following result is taken from [11], Theorem 3.8.ii. 

Proposition 1.4. Fix n > N , and choose parameters «j > for 1 < i < n and aj > for 
1 < i < N such that 7jj = — aj > 0. Consider the log-Gamma polymer partition function where 
dij ~ r _1 (7j j j). Then for all u with Re(u) > 



E 



e -uZ(n,N) 



s N (dw 1 , dw N ) ] [ T(aj - Wi) J 



.1 ' 



m N Coil f=i F w) 



where 

(1.2) F(w) = u w H T(a m - w). 



m=l 
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Until now, there has not been progress in extracting asymptotics from this formula. The fol- 
lowing theorem, however, transforms this integral formula into a Fredholm determinant for which 
we can readily perform asymptotic analysis. This identity should be considered the main technical 
contribution of this paper, of which Theorem 1 is essentially a corollary (after some asymptotic 
analysis) . 

Definition 1.5. We introduce the following contours: C$ is a positively oriented circle around the 
origin with radius 8; is is a line parallel to the imaginary axis from — ioo + 8 to ioo + 8; —is is 
similarly the contour from —ioo — 8 to ioo — 8; and i' Si 5a M is the horizontal line segment going from 
—8i + iM to 82 + iM. For any simple smooth contour 7 in C we will write L 2 {^/) to mean the space 
L 2 (7, /i) where \i is the path measure along 7 divided by 27ri. 

Theorem 2. Fix < 82 < 1, < 81 < min{<52, 1 — 82} and a\,... , ajv £ C such that \ai\ < 8\. 
Suppose F is a meromorphic function such that all its poles have real part strictly larger than 82, F 
is non-zero along and inside Cg 1 , and for all k > 
(1.3) 

/ dwe^-^ lm ^\lm(w)\ K \F(w)\<oo, [ dw e n ^- 1)llm(w)l \lm(w)\ K \F{w)\ ► 0. 

Then 

(1.4) / •••/ s N (d Wl ,...,dw N ) ' [T(a j -w l )]]-^- = det(I + K) L 2 { 



where det(J + K)^^ ) is the Fredholm determinant of K: L 2 (C<5 1 ) — > L 2 (Csj) with 



1 F(a,j) 

i,j=i j=i v ■>> 

2m, \ k t2i 



>- \ T-r, 1 f 7T F(w) 1 rr T(v — a m ) 

1.5 K(v,v') = — dw—- —y 

Z7ri J 1 sm(7r(u — it!) J -r (w) w — v 1 (it! — a m ) 

The proof of this theorem is given in Section 3. The argument uses only the Andreief identity 
[4] (sometimes referred to as the generalized Cauchy-Binet identity), the Cauchy determinant, the 
fact that det(I + AB) = det(I + BA) and some simple contour shifts and residue computations to 
go from the Fredholm determinant formula on the right hand side of (1.4) to the integral formula 
on the left hand side. Going from the integral formula to the Fredholm determinant formula (even 
knowing that it should be true) is a much more challenging path, since one has to undo certain 
cancelations to discover determinants. 

Remark 1.6. In the context of the related semi-discrete polymer (see the end of the introduction) 
O'Connell describes (after Corollary 4.2 of [14]) another Fredholm determinant expression for the 
Laplace transform of the partition function. That expression and the one which arises from the 
application of Theorem 2 are different. 

Remark 1.7. The condition in Theorem 2 on the \ai\ < 81 can be relaxed at the cost of more 
complicated choices of contours. Instead of taking the w contour to be a vertical line is 2 one could, 
in order to accommodate larger values of the a^, shift the vertical line horizontally to the right by 
some positive integer K, and then augment the w contour with a collection of sufficiently small 
contours around the positive integers up to and including K. See Section 5 of [8] for an example of 
this sort of procedure. 



We may apply Theorem 2 to show the following. 
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Corollary 1.8. Fix n > N and any 61,82 such that < 82 < 1 and < <5i < min{^2,l — 82}- 
Given a collection of on > 82 for 1 < i < n, and < a,- < 8\ for 1 < j < N, set 7$ = aj — a.,- 
and consider the log-Gamma polymer partition function with weights dij ~ r~ 1 (7j J ). Then for all 
u with Re(u) > 



E 



e -uZ(n,N) 



det(J + K u ) L 2 {Csi ), 



where det(I + K) L 2^ Cs ^ is the Fredholm determinant of K : L 2 (C<5 1 ) — > L 2 (C < 5 1 ) with 

r, ( A 1 / , 7T 1 -^r r(u-a m ) 

K u{v,v) = — dw—— vT^TT 7 J I TT7 7 

z7Ti sm(7r(z; — w)) I 1 (v) w — v' ±A - L [w — a m ) 



and F(w) as given in (1.2). 



Proof. We start with the Laplace transform formula given in Proposition 1.4 (if some a,j = we 
simply shift the iR integration contour slightly to the left). In order to apply Theorem 2 we must 
check that all the conditions are satisfied. The function F(w) has no poles with real part less 
than minjaj} and it is non-zero in the entire complex plane. This implies that given 8± and 82 
as specified in the hypothesis, the conditions on the poles and zeros of F are satisfied. It remains 
to check the decay condition (1.3). This, however, is immediate from the estimate on the Gamma 
function as given in (3.2). Note that the condition n > N becomes important in this case (actually 
n > N — 1 would do). The conditions having been satisfied, we may apply Theorem 2 to arrive at 
the corollary. □ 

The asymptotic analysis of this Fredholm determinant is performed in Section 2 and yields the 
proof of Theorem 1. 

The existence of such an identity as in Theorem 1 did not arise out of the blue. Let us briefly 
explain the two results which suggested this identity (though only in the special case of F(w) as in 
(1.6)). 

Definition 1.9. An up/right path in M. x Z is an increasing path which either proceeds to the right 
or jumps up by one unit. For each sequence < s± < ■ • ■ < sjv-i < t we can associate an up/right 
path eft from (0, 1) to (t, N) which jumps between the points (sj, i) and (sj, i + 1), for i = 1, . . . ,N— 1, 
and is continuous otherwise. Fix a real vector a = (a±,..., a^) and let B(s) = (Bi(s), . . . , B^(s)) 
for s > be independent standard Brownian motions such that Bi has drift Oj. 

Define the energy of a path (ft to be 

E((ft) = Bt^si) + (B 2 (s2) - B 3 (*l)) + • • • + (B N (t) - B N (8 N -i)) . 
Then the O'Connell-Yor semi-discrete directed polymer partition function Z N (t) is given by 

Z N (t) = [ d(fte E ^\ 



where the integral is with respect to Lebesgue measure on the Euclidean set of all up/right paths (ft 
(i.e., the simplex of jumping times < s\ < • • • < sat-i < t). 

Due to the invariance principle, the semi-discrete polymer is a universal scaling limit for discrete 
directed polymers when N is fixed and n goes to infinity (and temperature is suitably scaled - 
see for instance [5]). As such, the semi-discrete polymer inherits the solvability of the log-Gamma 
polymer. 
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In fact, before the work of [11] on the log-Gamma polymer, O'Connell [14] showed (see Corollary 



4.2) that E 



-uZ N (t) 



was given by the left hand side of (1.4) with 



(1.6) F(w) 



,.2 



t/2 



On the other hand, the theory of Macdonald processes was developed in [7]. Using Macdonald 
difference operators and the Cauchy identity for Macdonald polynomials [7] computes expectations 
of certain observables of the Macdonald processes, which, when formed into generating functions, 
lead to Predholm determinants. After a particular limit transition, the Macdonald processes be- 
come the Whittaker processes introduced in [14] and the generating functions become a Fredholm 



determinant expression for E 



-uZ N (t) 



(1.4). In particular the following resu 



- exactly in the form of the one on the right hand side of 
t appeared in [7] as Theorem 5.2.10 (a change of variables 



s + v = w brings it exactly to the form of (1-5)). 

Theorem 3. Fix N > 1 and a drift vector a = (a\, . . . , ajy)- Fix < 82 < 1, and 8± < 82/2 such 
that I ai| < 5\. Then for t > 0, 



-uZ N (t) 



det(J + K u ) 



E 

where det(I + K u ) is the Fredholm determinant of 

K u : L 2 (C a ) -> L 2 (C a ) 

for C a a positively oriented contour containing a±, . . . , a at and such that for all v, v' £ C a , we have 
\v — v'\ < 82 ■ The operator K u is defined in terms of its integral kernel 

1 r i r , , m , * rr r (" - <•»>) u-e""*" 2 ' 2 



, V + s 

771=1 

This provides a very indirect proof of the identity given in (1.4), in the very particular case of 
F(w) as in (1.6). An obvious question this development raised was to provide a direct proof of this 
identity and to understand how general it is. This is what Theorem 2 accomplishes. 

It is worth noting that Macdonald processes also have a limit transition to the Whittaker processes 
described in [11], which are connected to the log-Gamma polymer. The Fredholm determinant given 
in Corollary 1.8 was anticipated in [7] via the formalism of Macdonald processes yet its rigorous 
proof was so far lacking because of the nontriviality of certain decay estimates required by that 
approach. If one only cares about the GUE Tracy- Widom asymptotics then the approach given 
here provides a direct, though non-obvious and rather ad hoc, route from the integral formula of 
[11] to the Fredholm determinant. 

Acknowledgements. AB was partially supported by the NSF grant DMS-1056390. IC was par- 
tially supported by the NSF through DMS-1208998 as well as by the Clay Research Fellowship and 
by Microsoft Research through the Schramm Memorial Fellowship. DR was partially supported by 
the Natural Science and Engineering Research Council of Canada, by a Fields-Ontario Postdoctoral 
Fellowship and by Fondecyt Grant 1120309. DR is appreciative for MIT's hospitality during the 
visit in which this project was initiated. 

2. Fredholm determinant asymptotic analysis: Proof of Theorem 1 



Let us first recall two useful lemmas. In performing steepest descent analysis on Fredholm 
determinants, the following allows us to deform contours to descent curves. 
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Lemma 2.1 (Proposition 1 of [17]). Suppose s —>T S is a deformation of closed curves and a kernel 
L(rj,rf) is analytic in a neighborhood ofT s x T s C C 2 for each s. Then the Fredholm determinant 
of L acting on T s is independent of s. 

Lemma 2.2 (Lemma 4.1.38 of [7]). Consider a sequence of functions {f n }n>i mapping R — > [0, 1] 
such that for each n, f n (x) is strictly decreasing in x with a limit of ' 1 at x = — oo and at x = oo, 
and for each 5 > 0, on R \ [—5,5] /„ converges uniformly to l(x < 0). Define the r- shift of f n as 
fn{x) = f n (x — r). Consider a sequence of random variables X n such that for each r£l, 

nr n (x n )]^ P (r) 

and assume that p(r) is a continuous probability distribution function. Then X n converges weakly 
in distribution to a random variable X which is distributed according to ¥(X < r) = p(r). 

Proof of Theorem 1. Consider the function /at(x) = e e and define f^(x) = f n (x — r). Observe 
that this sequence of functions meets the criteria of Lemma 2.2. Setting 

u = n(n,r, 7 ) = e- n ^- rnl/3 

observe that 

~uz(n,n) = r ( log Z(n, n) -nf^ 

By Lemma 2.2, if for each rGlwe can prove that 

' log Z(n,n) - nL 



lim E 

n— »oo 



rr 
Jn 



nV3 

for p~f(r) a continuous probability distribution function, then it will follow that 

log Z(n,n) - n/ 7 



lim P — 7^ < r = p~(r) 

n— >oo V n i J 

as well. 

Now, the starting point of our asymptotics is the Fredholm determinant formula given in Corollary 
1.8 for E [e~ uZ ( n,n ^] . Observe that by setting aj = and Qj = 7, we can choose S± and 82 as necessary 
for the corollary to hold. In particular, that implies that 

'log Z(n,n) - n/ 7 s 



E 



rr 
Jn 



det(J + K u(jl)ra) ) 



nV3 

where det(/ + K u i n ^ r ^\) is the Fredholm determinant of K: L 2 (C < 5 1 ) — > L 2 (C < 5 1 ) with 

(2.1) K u{ ) ( v ,v') = -L[ dw—^ e^fniG^-Giw^+rn^iv-w))-^, 

v ' 2-iri j£ g sm(ir{v — w)) V / w — v' 

and 

G(z) = logr(z)-logr( 7 -z) + / 7 z. 

We have rewritten the kernel in the form necessary to perform steepest descent analysis. Let us 
first provide a critical point derivation of the asymptotics. We will then provide the rigorous proof. 
Besides standard issues of estimation, we must be careful about manipulating contours due to a 
mine-field of poles (this ultimately leads to our present technical limitation that 7 < 7*). 

The idea of steepest descent is to find critical points for the argument of the function in the 
exponential, and then to deform contours so as to go close to the critical point. The contours should 
be engineered such that away from the critical point, the real part of the function in the exponential 
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decays and hence as n gets large, has negligible contribution. This then justifies localizing and 
rescaling the integration around the critical point. The order of the first non-zero derivative (here 
third order) determines the rescaling in n (here n 1//3 ) which in turn corresponds with the scale of 
the fluctuations in the problem we are solving. It is exactly this third order nature that accounts 
for the emergence of Airy functions and hence the GUE Tracy- Widom distribution. 

Let us record the first three derivatives of G(z): 

G'{z) = *(z) + *( 7 -z) + / 7 , 

G"(z) = m'(z)-^'{ 1 -z), 

G"'{z) = V"(z) - *"(7 - z). 

Define t 7 such that G"(t~) = 0; that is, t~ = 7/2. Notice that f~ f was defined such that G'(t~ f ) = 
as well. Around z = t 7 , the first non- vanishing derivative of G is the third derivative. Notice that 
g 7 = —G"'(tJ). This indicates that near z = f~, 

G(v) - G(w) = _ frfr-*r) 3 + fr(" -*r) 8 + h . . t ., 

6 6 

where h.o.t. denotes higher order terms in (v — t 7 ). This cubic behavior suggests rescaling around 
t 7 by the change of variables 

v = n 1/3 (v - w = n l/3 (w -t 7 ). 

Clearly the steepest descent contour for v from i 7 departs at an angle of ±2-7r/3 whereas the contour 
for w departs at angle ±7r/3. The w contour must lie to the right of the v contour (so as to avoid 
the pole from l/(w — v')). As n goes to infinity, neglecting the contribution away from the critical 
point, the point-wise limit of the kernel becomes 

, /N 1 f 1 exp{^v 3 +rv\ M 
(2.5) K rn {v,v')- 1 



where the kernel acts on the contour e~ 2n1 ' 3 M+ U e 27r1 ' 3 M+ (oriented from negative imaginary part 
to positive imaginary part) and the integral in w is on the (likewise oriented) contour {e -7rl / 3 K+ + 
6} U {e 7ri ' 3 R4. + 6 for any horizontal shift 5 > 0. This is owing to the fact that 

1 /o 7r 1 dw dw 

n ' > > . 

sin(7r(i> — w)) v — u) ' w — v' w — v' 

where the n -1 ' 3 came from the Jacobian associated with the change of variables in v and v' . 
Another change of variables to rescale by (^ 7 /2) 1 / 3 results in 

K (vv')= 1 / 1 ex P(-^ 3 / 3 + (g7/ 2 )" 1/3 ^} dw 
r ' 7 ' 27ri J v — w exp {— w 3 /3 + (g 1 /2)~ 1 / 3 rw] w — v' 

One now recognizes that the Fredholm determinant of this kernel is one way to define the GUE Tracy- 
Widom distribution (see for instance Lemma 8.6 in [8]). This shows that the limiting expectation 
p 7 (r) = -Fgue ((57 /2)~ 1 / 3 r) which shows that it is a continuous probability distribution function 
and thus Lemma 2.2 applies. This completes the critical point derivation. 

The challenge now is to rigorously prove that det(I + K u i n ^^\) converges to det(I + K ra ) where 
the two operators act on their respective L 2 spaces and are defined with respect to the kernels above 
in (2.1) and (2.5). 
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We will consider the case where 7 < 7* for 7* small enough and perform certain estimates given 
that assumption. Let us record some useful estimates: for z close to zero we have {g is being used 
to represent the Euler Mascheroni constant g = 0.577...) 

logr(z) = — log z — gz + 0(z 2 ) 

= -l-g + 0(z) 

= \ + o{i) 

z z 

= -^ + 0(1) 



6 



= f 4 +0(l). 
From this we can immediately observe that for 7 small, 

h = 7/2 

A = 4 7 - 1 + 2 5 + 0( 7 ) 
g 7 = 32 7 - 3 + 0(l). 

With the change of variables z = 7Z we may estimate 

G(z) - G(i 7 ) = logics) - lo g r( 7 - z) + / 7 (z - 7 /2) + - f K t K 
= f(z) + 0( 1 2 ), 
where the error is uniform for z in any compact domain and 

f{z) = log(l - z) - log(z) + 4z - 2. 

Our approach will be as follows: Step 1: We will deform the contour Cs 1 on which v and v' 
are integrated as well as the contour on which w is integrated so that they both locally follow the 
steepest descent curve for G{z) coming from the critical point i~ and so that along them there is 
sufficient global decay to ensure that the integral localizes to the vicinity of i~. Step 2: In order 
to show this desired localization we will use a Taylor series with remainder estimate in a small ball 
of radius approximately 7 around t~, and outside that ball we will use the estimate (2.9) for the 
v and v' contour, and a similarly straightforward estimate for the w contour. Step 3: Given these 
estimates we can show convergence of the Fredholm determinants as desired. 

Step 1: Define the contour Cf which corresponds to a steep descent contour 1 of the function f(z) 
leaving z = 1/2 at an angle of 27r/3 and returning to z = 1/2 at an angle of — 27r/3, given a 
positive orientation. In particular we will take Cf to be composed of a line segment from z = 1/2 to 
z = 1/2 + e 27rl//3 , then a circular arc (centered at 0) going counter-clockwise until z = 1/2 + e~ 2?n / 3 
and then a line segment back to z = 1/2. It is elementary to confirm that along Cf, Re(/) achieves 
its maximal value of at the unique point 2 = 1/2 and is negative everywhere else. The contour 
7C7 will serve as our steep descent contour for G(v) (see Figure 1). 

One should note that for any r G (0, 1), for 7 small enough, the contour jCf is such that for v 
and v' along it, \v — v'\ < r. By virtue of this fact, we may employ Proposition 2.1 to deform our 
initial contour C$ ± to the contour ^Cf without changing the value of the Fredholm determinant. 



^For an integral I = J dze ti( - z \ we say that 7 is a steep descent path if (1) Re(/(z)) reaches the maximum at 
some Zo £ 7: Re(/(z)) < Re(/(zo)) for z£) \ {zo}, and (2) Ke(f(z)) is monotone along 7 except at its maximum 
point zo and, if 7 is closed, at a point z\ where the minimum of Re(/) is reached. 



LOG-GAMMA POLYMER FREE ENERGY FLUCTUATIONS VIA A FREDHOLM DETERMINANT IDENTITY 9 




Figure 1. Steep descent contours 

Likewise, we must deform the contour along which the w integration is performed. Given that 
v,v' £ 7C7 now, we may again use Proposition 2.1 to deform the w integration to a contour C/ n 
which is defined symmetrically over the real axis by a ray from i 7 + n -1 / 3 leaving at an angle of 7r/3 
(again, see Figure 1). It is easy to see that this deformation does not pass any poles, and with ease 
one confirms that due to the linear term in G(w) the deformation of the infinite contour is justified 
by suitable decay bounds near infinity between the original and final contours. 

Thus, the outcome of this first step is that the v,v' contour is now given by 7C/ and the w 
contour is now given by C/ n which is independent of v. 

Step 2: We will presently provide two types of estimates for our function G along the specified 
contours: those valid in a small ball around the critical point and those valid outside the small ball. 
Let us first focus on the small ball estimate. 

Lemma 2.3. There exists 7* > such that for all 7 < 7* the following two facts hold: 

(1) There exists a constant c\ > such that for all v along the straight line segments ofjCf: 

Re [G(v) - G(i 7 )] < Re [-a-f 3 {v - t 7 ) 3 ] . 

(2) There exists a constant C2 > such that for all w along the contour C/ n at distance less than 
7 from t 7 : 

Re [G(w) - G(ty)] > Re [-c 2 ^{w - t 7 ) 3 ] . 
Proof. Recall the Taylor expansion remainder estimate for a function F{z) expanded around z, 
\F(z) - (F(z) + F'(z){z -z) + \F"{z)(z - zf + \F'"{z)(z - zf) \ < max ^\F'"> {Q\\z-z\\ 

Q£B(z,\z-z\) 



LOG-GAMMA POLYMER FREE ENERGY FLUCTUATIONS VIA A FREDHOLM DETERMINANT IDENTITY 10 



We may apply this to our function G(z) around the point t 7 giving 

\G(z) - Gfo) + l^(z - i 7 ) 3 | < max _ ±\G""(()\\z - t^. 

C&B(ty,\Z— ty\) 

Let z = jz and also let £ = 7^. Then for 7 small, we have the following estimate 

\G(z)-G(t,) + f(z-l/2f\<\ 

Both parts of the lemma follow readily from this estimate and the comparison along the contours 
of interest of -§-(z — 1/2) 3 with the right hand side above. □ 

We may now turn to the estimate outside the ball of size 7. 

Lemma 2.4. There exists 7* > and c > such that for all v along the circular part of the contour 
^Cf, the following holds: 

Re[G(v) - G(ly)] < -c. 

Proof. Writing v = "fv we may appeal to the second line of the estimate (2.9) and the fact that 
along the circular part of the contour "yCf, the real part of function f(v) is strictly negative and 
the error in the estimate is 0(7). □ 

The above bound suffices for the v contour since it is finite length. However, the w contour is 
infinite so our estimate must be sufficient to ensure that the contour's tails do not play a role. 

Lemma 2.5. There exists 7* > and c > such that for all w along the contour C/ n at distance 
exceeding 7 from i 7 , the following holds: 

Re[G(iu) - G(tj)] > Re[c7 -1 u;] . 

Proof. This estimate is best established in three parts. We first estimate for £ between distance 7 
and distance C17 from f 7 (c\ large). Second we estimate for w between distance C17 and distance C2 
from t^,. Finally we estimate for all w yet further. This third estimate is immediate from the first 
line of (2.9) in which the linear term in z clearly exceeds the other terms for 7 small enough and 
\z\ > c 2 . 

To make the first estimate we use the bottom line of (2.9). The function f(z) has sufficient 
growth along this contour to overwhelm the 0(7 2 ) error as long as 7 is small enough. The 0(7 2 ) 
error in (2.9) is only valid for z in a compact domain though. So, for the second estimate we must 
use the cruder bound that G(z) — G(i 7 ) = f(z) + 0(1) for z along the contour and of size less than 
C2 from t^. Since in this regime of z, f(z) behaves like (z — 7/2)/ 7 = A'y~ 1 z + 0(1), one sees that 
for 7 small enough, the O(l) error is overwhelmed and the claimed estimate follows. □ 

Step 3: We now employ the estimates given above to conclude that det(I + K u i ntT ^yj converges to 
det(/ + K r ^ry) where the two operators act on their respective L 2 spaces and are defined with respect 
to the kernels above in (2.1) and (2.5). The approach is standard (see for instance [8, 3]) so we 
just briefly review what is done. Convergence can either be shown at the level of Fredholm series 
expansion or trace-class convergence of the operators. Focusing on the Fredholm series expansion, 
the estimates provided by Lemmas 2.3, 2.4 and 2.5, along with Hadamard's bound, show that for 
any e there is a A; large enough such that for all n, the contribution of the terms of the Fredholm 
series expansion past index k can be bounded by e. This localizes the problem of asymptotics 
to a finite number of integrals involving the kernels. The same estimates then show that these 
integrals can be localized in a large window of size n" 1 / 3 around the critical point t 7 . The cost of 



z 3 



(i-~zY 



l/2| 4 + 0( 7 2 ). 
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throwing away the portion of the integrals outside this window can be uniformly (in n) bounded by 
e, assuming the window is large enough. Finally, after a change of variables to rescale this window 
by n 1 / 3 we can use the Taylor series with remainder to show that as n goes to infinity, the integrals 
coming from the kernel K u ^ n ^ r ^ converge to those coming from K r>T This last step is essentially 
the content of the critical point computation given earlier. □ 



3. Equivalence of contour integral and Fredholm determinants 



Let us first recall a version of the Andreief identity [4] (sometimes referred to as the generalized 
Cauchy-Binet identity): Assuming all integrals exist, 



(3.1) 



1 

aT 



dfi(zi) • • • dfi(z N ) det [fi(zj)]^ j=1 det [gi(zj)]^ j= 



det 



dfJ>(z) fi(z)gj(z) 



N 



Proof of Theorem 2. We start working from the right hand side of the formula. We regard the 
Fredhom determinant as given by its Fredholm series expansion (the convergence of the series is 
a consequence of the identity we will prove). Note that our assumptions on 8\ and 82 imply that 
8\ < |, so the factors T(v — a m ) have no poles for v G Cg 1 . They also imply that the factor 
sin(7r(-y — to)) (to — v') in the denominator in (1.5) is never zero on £$ 2 . To see that the integral is 
convergent, recall the asymptotics \T(x + iy)\e n ^ 2 ^ |y| 1//2_x — > \p2m as y — > ±00 for x,y G R (see 
(6.1.45) in [1]), which implies that, for Re(z) in some finite interval, 



(3.2) 
as |Im(z 



I — > 00 for some c\ > and i] G 
sml7r,2)| > C2e 7r ' lm ^)' for some C2 > 0. Then 



T(z)\ > cie-fl Im WI|Im(z)| ? ' 

We also have, under the same assumption on z, 



(3.3) 



7T 



F(w) 1 



sin(-7r(t> — w)) F(v) w — v ' 



N 

n 

m=l 



T(w - a m ) 



< 



c 3 e 



^(f-l)|Im(«,)|| Im ^|-7V^-l| F ^| 



as |Im(ti))| — > 00 for some C3 > 0, and hence the integral in (1.5) is convergent by (1.3). 
Observe that K = AB with A: L 2 (4 2 ) — > L 2 (C Sl ) and B:L 2 (C Sl ) — > L 2 {£ &2 ) given by 



A(v, to) 



7T 



N 



sin(-7r(t; — to)) F(v) 



F(w) yj T(v - a m ) 
(w - a m ) 



m=l 
2 



and 



B(w,v') 



to 



(checking that A and B map L (£g 2 ) to L 2 (C$ 1 ) and vice- versa is easy, and uses (1.3) in the case of 



A). Let K = BA: L 2 (l 



S2 ■ 



L (£$ 2 ), which is then given by 



K(w, to') = — : 



2vri 



dv 



7T 



F(w') 



Csi w-v sin(7r(t> - w')) F(v) T(w' - a m ) 



N 

n 



T(v 



Now recall in general that if A and B have integral kernels so that the integrals j dw A(v, w)B(w, v') 
and J dv B(w,v)A(v,w') both converge for v,v',w,w' in the right domains, then det(J + AB) = 
det(I + BA) in the sense that the two convergent Fredholm expansions coincide termwise (this can 
be seen as an application of the Andreief identity (3.1)). Using this we deduce that 



(3.4) 



det(I + K) L 2 



det (J + K) L 2 



(tsa)- 
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N 



Recalling that T(l + s) = sT(s), let 
(3.5) G(v) = F(v) H f 

and rewrite K as 



v — a r . 



1 T(v - a m + 1) 



K(w, w') = — — 



dv 



1 



7T 



G(w>) 



2vri Jc s w — v sin(7r(t; — it/)) G{v) 

By the assumption 5\ < min{<52, 1 — £2}, w — v and sm(ir(v — w')) are never zero for v £ Cs 1 and 
10,10' 6 £^ 2 . Thus the only poles of the integrand inside C$ 1 are those coming from G(t>) _1 , which 
can only come from the factors (v — a m ) in the definition of G(v). For simplicity we will assume 
that all the a^s are pairwise distinct, the general case follows likewise by taking limits. Then the 
integrand has simple poles at v = cii, i = 1, . . . , N, and thus 



K(w, w') = —ir G(w')Res 

i=l 

1 N 1 



7T 



w — v sm(ir(v — w')) G(v) 

7T _L -i-r T(di - a m + 1) 



i=l 



We rewrite this as 
with 

(3.6) fi(w) 



w — cii sin(7r(aj — w')) F(cii) 



1 N 

K(w,w') = —^2fi{ w )9i(w' 



a; - a r 



w — Cli 



9i{w 



CiG{w') 



i=l 



7T 



sin(7r(aj — w')) 



and Cj 



n r ^ - a * 



In particular, this means that K = AB where A : £ 2 ({1, ■ ■ ■ , N}) — > L 2 (£g 2 ) is given by A(w, i) = 
(2m)- 1 fi(w) and B: L 2 {£ S2 ) ({I,..., N}) by B(i, w') = 9i (w'), so that 

BA(i,j) = —f dw fi(w)gi(w). 

The integral is convergent for the same reason as the integral in (1.5), and hence using again the 
formula det(7 + AB) = det(7 + BA) we deduce that 

r i N 

det(I + K) L 2 (ig2) =det 



1 



li=j + — dwfi(w)gj(w) 
2m Jis 2 



»J=1 



Our next goal is to shift the contour £$ 2 in the integral appearing on the right hand side of (3.7) 
to —£$ 1 . Note that, as we do this, we will cross all the etj's. We have 



fi(w)gj(w) = F(w) 



CjTT 



sin(7r(aj — w)) 



1 n 1 

r(to - Oi + 1) A J-. T(w - 



a-m) 



Since all the poles of F lie to the right of £g 2 , the only singularities of fi(w)gj(w) we encounter as 
we shift the contour are those coming from the zeros of the sine, and hence we only see simple poles 
at w = en in the case i = j. On the other hand, the integral of fi(w)gj(w) along the segments going 
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from —5i ± \M to 62 ± iM goes to as \M\ — > 00 by (1.3) and (3.3) as before. The conclusion is 
that 

— / dw fi(w)gj(w) = l i=j Res [fi(w)gi(w)} + — [ dw fi(w)gj(w). 

27T1 J e W=CH 27T1 J_ £ 



Since 



Res[f i (w)g i (w)] = -C i F(a i )\ 



a,- - a r . 



W = CLi 



T(ai - a m + 1) 



1. 



we deduce from (3.7) that 

det(I + K) L 2 {£h) = det 



(3.8) 



1 
2^i 



dw fi(w)gj(w) 



N 



J i,j = l 



(27ri) ljV iV! J e I e dwi "' dwN det iM w j)]Z=i de %K)]i 



N 

j=l, 



where the second equality follows from the Andreief identity (3.1). 

Note that the contours in this last integral are exactly the ones appearing in the formula we 
seek to prove. Hence what is left is to compute the two determinants and show that the integrand 
coincides with the one appearing in (1.4). We start by observing that det[fi(wj)]fj =1 is a Cauchy 
determinant, so 



(3.9) 

On the other hand 



det[/i(uy)]£. =1 



n»<j( a i - ai)(wj - wj] 



N 



det[5i(u;j)]ij = i = Y[G(wi)Cidet 



i=l 



7T 



sin(7r(ai — Wj)) 



N 



The determinant on the right hand side can be turned into another Cauchy determinant by writing 

7r 27ri 1 



Sin(7r(aj — Wj)) e -i7r(ai+Wj) e 



m{ai+Wj) p2mai „2mwj 



so that 



det 



7T 



sin(7r(aj — Wj)) 



N N . 

nZTTl 
— ^— — rdet 



i,j=l i=l 



-iiz(ai+Wi) 



1 



g2i7rai ^2mujj 



N 



and therefore 

(3.10) det^(^)]£=i = (2^)^^+^ J] G(wi)Ci 



N 



i=l 



J-J ^ e 2inaj _ e 2ina,i^ e 2inwi _ e 2mwf 

Y[ N 1 (e 2i7rai — e* mw i~) 



Now we collect the factors involving only Wi's, only aj's, and cross-terms: from (3.9), (3.10) and the 
definitions (3.5) and (3.6) of G(w{) and Cj we get 



(3.11) 



detlfiiw^j^detlgiiwj)}^ = (2iri) N D a , a D W:W D a 
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with 

D a , a = e^- H( aj - a t )(e 2 ™> - e 2 ^) {[ ""^^ " ^ ' 

i<j i=l 
N 

D w , w = e i7r ^ Y[F( Wi ) Y[( Wi - w 3 ){e 2 ^ - e 2i ™>), 

i=l i<j 

d — TT - 1 1 

a,w J- -I- Wi - a,- T(wi - aA e 2i7ra * - e 2inw j ' 

i,j = l J 

a can be rearranged in the following way: 

D aa = e in ^ a * TT(aj - a i )e i7r(0i+a ' ) (e i7r( ° i ~ a, ' ) _ e in( aj - ai ) ^i^r V ( Ql ~ a ^ 

AT 

= e-JVE, ^ ( 2i )|^-D Jj r(a . _ + i)r( a . _ ai ) s i n (^( ai _ q.)) JjF(oi)- 1 

i<j i=l 
iV 

= e ^E i a i (_2 7ri )|iV(Ar-i)-Q j p (a . ) -i ) 

i=l 

where have used 

r(_ s )r(i + s ) 



sin(— 7rs) 
Similar computations give 

A.,* = e^S^ l(27ri) ^(^-i) jj M<m- Wj )) {wi _ 

i<j i=l 



N 

D a , w = e-™ N ^+^)(27Ti)- N2 J] r(a 

Hence 



D a ,a,D w w D a ,w 



1 , ^ijv(jy-i) tt sin(7r(w, - w 3 -)) , y-r . -p-j- F(wj) 

v ' i<j ij=l i=l v v 

1 1 □ 



(27Ti)^iir(^-^)l = lF(G l )' 

and thus from (3.8) and (3.11) we deduce that 



~ f f N N p^ w ,\ 

det(/ + K) LHl6i } = y _^ • • • dm- dw N J nw ._ Wj) U T ^ -m)Up^- 

Putting this together with (3.4) and the definition (1.1) of sat yields (1.4). □ 
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